
load datos3a.mat
x = log(abs(xd));
y = log(abs(yd));

P = polyfit(x,y,1);
m = P(2);
b = P(1);

A = m;
C = e^b;

plot(xd,yd,'r',xd,C*(xd.^A),'b')

